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Abstract 

The three-dimensional nonlinear dynamics of an electron gas in a 
semiconductor quantum well is analyzed in terms of a self-consistent 
fluid formulation and a variational approach. Assuming a time-de- 
pendent localized profile for the fluid density and appropriated spatial 
dependences of the scalar potential and fluid velocity, a set of ordinary 
differential equations is derived. In the radially symmetric case, the 
prominent features of the associated breathing mode are characterized. 

1 Introduction 

The fast dynamics of a quantum electron gas confined in nanoscale systems 
including quantum wells, quantum dots and metal clusters is receiving a great 
deal of interest [U [2] . Among possible collective motions, one has the electric 
dipole response, the so-called Kohn mode [3], associated to rigid oscillations 
of the electron gas at the trapping frequency. 

Other fundamental collective oscillations are the breathing modes, char- 
acterized by a pulsating electron cloud. Such breathing modes have been 
investigated through numerical simulation of the N-body time-dependent 
Schrodinger equation for restricted number of particles, exploring large en- 
ergy coupling parameter regimes [3]. Soon later, a variational approach based 
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on a mean-field, quantum hydrodynamic model has been used to analyti- 
cally describe the nonlinear breather dynamics of an one-dimensional electron 
gas in a semiconductor quantum well [5J. A time-dependent Rayleigh-Ritz 
method was employed in order to reduce the underlying quantum fluid model 
to a set of nonlinear second-order ordinary differential equations. The nu- 
merical analysis of the kinetic, Wigner-Poisson system was shown to be in 
agreement with the analytical results. In the present work we make a first 
step toward the extension of this approach to three spatial dimensions. 



2 Model equations 

In the mean-field approximation, the electron dynamics can be described by a 
self-consistent quantum hydrodynamic model that was also applied for quan- 
tum plasmas [6 J and metallic nanostructures [7] . Here only electrostatic fields 
are considered. Then the evolution of the electron fluid density n(r, t) and 
the electron fluid velocity u(r, t) follows from the continuity and momentum 
equations, 

^ + V-(nu) = 0, (1) 

— + u-Vu = + ^^ V — 7=- » 2 

at m*n m* zm, \ yn J 

where m* is the effective electron mass, h Planck's constant divided by 2tt, 
P = P(r, t) the electron fluid pressure and 

V eS = V coni (r) + V H (r,t) (3) 

the effective potential, which is composed by a confining V con f and a Hartree 
Vh terms. The Hartree potential obeys Poisson's equation, 

W fl = , ( 4 ) 

where e is the magnitude of the electron charge and e is the effective dielectric 
permeability in the semiconductor quantum well. The term proportional to 
h 2 on the right-hand side of Eq. (j2J) is associated to the so-called Bohm 
potential and is responsible for the quantum wave-like effects. 

For the confinement suppose the generally anisotropic potential 

2 

Vconf = „ {KlX + K 2 y + K 3 Z ) , (5) 
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where ^1,2,3 > are dimensionless constants, not necessarily equal. Moreover, 
the frequency uo can be related to a fictitious homogeneous positive charge 
of density n through co = (e 2 no/ m^e) 1 ' 2 . 

It is convenient to normalize time to coq 1 , space to Lq = (ks T /m*) 1 ' 2 /ujq 
[where ks is Boltzmann's constant and T the temperature], density to rio, 
velocity to L u>o, pressure to n fee T and energy to ks T. These are the same 
rescaled variables as in Refs. [HE)]. For simplicity, old and new variables will 
be represented by the same symbols. The model equations then become 

^ + V- (nu) = 0, (6) 

„ VP ^ „ Tr H 2 ( V 2 y/n\ 
— + u-Vu = Kt-VVh + — V — ^- , (7) 



dt n 

W H = -n, (8) 

where K is the diagonal dyad with components = S^Kj. Quantum effects 
are measured by the dimensionless parameter 

In addition, in the new coordinates we have 

Konf = HlX 2 + K 2 y 2 + K 3 Z 2 . (10) 

In order to close the system, the pressure in Eq. (J7|) must be related to 
the electron fluid density. In the fast time-scale implied by the trapping fre- 
quency, an adiabatic equation of state is indicated, since there is no sufficient 
time for thermalization. In this regard, the temperature parameter T can be 
interpreted as a measure of the average kinetic energy per electron. Hence 
we choose the polytropic relation 
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(11) 



where 7 = 5/3 is the three-dimensional polytropic exponent, and n a refer- 
ence density, whose choice will be discussed later. 

Let us make an estimate of the parameters in the model. For semi- 
conductor quantum wells [8], typically we have: the effective electron mass 
and the effective dielectric permeability are, respectively, m* = 0.067m e and 
e = 13 Eq and the equilibrium density is n = 4.7 x 10 22 m -3 . These values 
yield an effective plasmon energy huo = 8.62 meV, a characteristic length 
Lq = 16.2 nm, a Fermi temperature Tp = 51.8 K, and a typical time scale 
cuq 1 = 76 fs. An electron temperature T = 200 K then corresponds to a value 
if =0.5. 
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3 Time-dependent Rayleigh-Ritz method 



In order to derive a closed system of ordinary differential equations in which 
time is the independent variable, we first express the quantum hydrodynam- 
ical equations in a Lagrangian formalism. This approach is not based on a 
perturbative expansion, and is thus not restricted to the linear regime. A 
Lagrangian density corresponding to the system read as 



1 d6 r n 

£ = ^V H f -nV eff - n—-J W(n')dn' 



(12) 



where the independent fields are n, 9, and Vh- The velocity field follows from 
the auxiliary function 9 = 9(r,t) through 



u 



V0. 



The quantity W(n) in Eq. ffl2l) originates from the pressure, 



W = 



n dP dn' 
dn' n' 



(5/2) (n/n 



■rrW3 



(13) 
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Taking the variational derivatives of the action S = J Cdrdt with respect 
to n, 9, and Vh, we obtain the Eqs. (E|) (J8]) . 

The existence of a variational formalism can be used to derive approxi- 
mate solutions via the time-dependent Rayleigh-Ritz trial-function method. 
The electron density can be taken in the form of an anisotropic Gaussian, 
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while the self-consistent electrostatic field can be taken as 

Erf [R/y/2) 



(15) 



V, 



H 



,1/3 



R 



(16) 



where the variable 



R 



x — d\ 



y - d 2 



1/2 



+ 
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(17) 
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was introduced. Here <7j and di are time-dependent functions respectively 
related to the width and center-of-mass coordinate of the electronic cloud. 
In addition, A is a numerical constant directly related to the number of 
electrons in the quantum well, 

Finally, Erf is the error function, defined by 




The arguments in favour of the selected profile are: it describes a local- 
ized solution amenable to relatively simple analytical calculations; the ground 
state for the three-dimensional anisotropic quantum harmonic oscillator is a 
Gaussian. Hence we have an exact solution for negligible Hartree and pres- 
sure terms; in the radially symmetric case (<7i = o<i = 03) the proposed form 
is an exact solution for Poisson's equation. According to the last reasoning, 
in principle the model becomes less accurate in strongly anisotropic cases. 

To complete the Ansatz, we need to specify the velocity field, which we 
choose so as to exactly solve the continuity equation for the given n. Hence 
we take the linear form 

«i = -(r i -d i ) + d i , i = 1,2,3, (20) 

for the i— component of the velocity field. 

Since u = W, the variable 6 in the Lagrangian density can be written as 

e = E <r* ~ rf ^ 2 + * ( r * - *)) ( 21 ) 

ignoring an irrelevant additive purely time-dependent contribution. 



(18) 



4 Mechanical system 

Now a system of Newton equations for a iy di can be derived. Apart from a 
multiplicative factor, the Lagrangian turns out to be 

L = J CdY= l -a 2 + ^d 2 -U a (a)-U d (d), (22) 
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where a = (o"i, cr 2 , <7 3 ), d = d 2 , d 3 ) and 
U a = U a (a) = -(niaf + K 2 crl + K 3 a 2 ) 

A(a 1 a 2 a 3 ) 1 ' 3 ( 1 1 l\ 9 [s ( A\ 2/3 1 
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U d = U d (d) = ^(K 1 d 2 1 + K 2 dl + K 3 d 2 3 ) (24) 

are pseudo-potentials resp. for the breathing and dipole motions. 

Examining JJ a we find: the term with Ki )2 ,3 is due to the harmonic con- 
finement; the ~ H 2 term is due to the Bohm potential; the A terms are 
due to the self-consistent Hartree energy; the ~ (A/n) 2 ^ 3 term is due to the 
thermodynamic pressure. 

Clearly the dipole and electronic cloud width dynamics are decoupled, a 
feature of harmonic traps. The dipole motion is linear, 

d] + Kidi = , i = 1, 2, 3 , (25) 

corresponding to Kohn oscillations. On the other hand, the Oi execute cou- 
pled nonlinear oscillations. 

The equations for the breathing motion read 

V2^4_ AS A S 3 1 

ai + Ki<Ti ~ ^J^ l ~3V27 3 + W2^h^ 



H 2 /3\ 3 / 2 /A\ 2/3 1 



where 



S = S(t) = (<x 1( x 2 a 3 ) 1/3 . (27) 

A more detailed investigation of the dynamical system f|26|) will be post- 
poned to future work. Instead, in the following the radial case is considered. 

5 The isotropic case 



Assuming ^1,2,3 = l,cr 12) 3 = er, the breather equation reduces to 

A (H 2 /3\ 3 / 2 (A\ 2 ' 3 \ 1 



Before proceeding, we can now define the parameter n in the equation of 
motion, and we do this using thermodynamic arguments. In terms of our 
Gaussian Ansatz, we can expect 

n ~ -3 , (29) 

where o"o is the equilibrium value of the variance a. Hence A/n is not negli- 
gible even if A — > 0, which means neglecting the self-consistent interaction. 
Taking the simultaneous limits A — > 0, H — > with fixed A/n we have 



/ 3\ 3 / 2 /A\ 2/3 1 



On the other hand, restoring physical coordinates and assuming energy 
equipartition, we have the ensemble averages 

m^uol < x 2 > _ m^u 2 , < y 2 > _ m^tu 2 < z 2 > _ k b T 

2 ~ 2 ~~ 2 ~^T' ^ ' 

which can be also thought as the definition of the temperature parameter T. 
Supposing 

< x 2 > = < y 2 > = < z 2 > = al (32) 
we find cr = L . Going back to the rescaled variables, we expect 

cto = 1 (33) 

to be the equilibrium value in accordance with energy equipartition, neglect- 
ing Hartree and quantum terms. On these grounds, from Eq. (1301 we have 



/3\ 9/4 

n = ( - J A ~ 0.32A (34) 

In a first approximation, the above choice is adopted also for the A ^ 
0, H ^ case. Therefore we get 

A ( H 2 \ 1 
6V2a 2 \ J 

or 

dU 



Figure 1: Qualitative shape of the pseudo-potential U(a) for arbitrary A > 
0,H>0 



where U(a) is a pseudo-potential defined by 



a 2 



From the shape of the pseudo-potential (see Fig. [Q, it follows that a 
will always execute nonlinear oscillations around the unique minimum cr = 
<To(A,H), which is a solution of the algebraic equation U'(ao) = 0, or 

1-— + — ) 4 = 0. (38) 

The frequency Q = Q(A, H) of the breather mode, corresponding to the 
oscillations of a, can be obtained by linearizing the equation of motion in the 
vicinity of the minimum of U(a). In other words, we use o"o from Eq. (1381) 
plugging it in 

„ (d 2 U\ A ( H 2 \ 1 

A 

= 4- 



the last equality following from Eq. fl38|) and showing that — > 2 in the 
absence of self-consistent interaction. 
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A consistency check can be performed neglecting the Hartree potential, 
which yields the equilibrium variance 




a ° = I 1 + T ) ■ (40) 
This expression displays the correct low- and high-temperature limits for 



the quantum harmonic oscillator: cr — > 1 for H — > 0; and cr ~ yH/2 for 
H » 1. These estimates can be found from standard calculations on the 
canonical ensemble of quantum harmonic oscillators 

Since from Eq. ( |39|) we automatically have Q 2 > we conclude that 

°l > —r . ( 41 ) 

which further confirms the repulsive effect of the Hartree potential. 

On the other hand, Eq. (|38|) shows that for large Hartree parameter we 
have 

ff «~671- (42) 

Using this in Eq. ( 1391) we find that £1 ~ \/3 for large A. 
We can define 

f n 2 dr 1 A 
<n>= J —— = —= 43 
/ nc/r 2^2 erg 

as the average electron fluid density, where the average is calculated using 
the density iteself, at the equilibrium value a = ctq. Using n for calculating 
expectation values is reasonable since it corresponds to the quantum proba- 
bility density function. 

Table [1] shows that the breather frequency depends weakly on the param- 
eter H (and hence on the electron fluid temperature). The density < n > 
is shown to decrease with H. This follows since quantum diffraction effects 
tend to enlarge the width of the electronic cloud, which for fixed number of 
electrons means a smaller density. 



Table [2] shows the dependence of the breather frequency on the Hartree 
parameter A, for fixed H = 0.5 which is representative of realistic semicon- 
ductor quantum wells. In addition, the average particle density is shown. 
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H 


00 


n 


< n > 


0.0 


1.03 


1.97 


0.32 


0.5 


1.04 


1.97 


0.31 


1.0 


1.08 


1.98 


0.28 


1.5 


1.14 


1.98 


0.24 


2.0 


1.21 


1.98 


0.20 


2.5 


1.28 


1.99 


0.17 


3.0 


1.36 


1.99 


0.14 



Table 1: Equilibrium width ctq, breather frequency f2 and mean particle 
density < n > for A = 1 and various values of H 



A 


0"o 


n 


< n > 


0.0 


1.02 


2.00 


0.00 


1.0 


1.04 


1.97 


0.31 


2.0 


1.07 


1.95 


0.58 


3.0 


1.10 


1.93 


0.80 


4.0 


1.12 


1.92 


1.00 



Table 2: Equilibrium width a , breather frequency f2 and mean particle 
density < n > for H = 0.5 and various values of A 

6 Conclusion 

In this work, we derived the basic equations for the three-dimensional varia- 
tional description for the many-electron of a quantum electron gas in a semi- 
conductor quantum well. The first results on the associated linear breathing 
frequencies were then obtained. The analytic calculations should be con- 
fronted against numerical simulation of the Wigner-Poisson system, as made 
in the one spatial dimension case in [5]. In particular, the adequacy of the 
Gaussian profile proposed here in the case of large Hartree energy should be 
investigated. Moreover, the effect of the coupled dipole-breather dynamics in 
non-parabolic wells is also an interesting issue, since in this case the breather 
mode can be triggered using a purely dipolar excitation [5]. Standard pump- 
probe experiments can then be used to optically detect the breather mode. 
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